Specific heat behavior of high temperature superconductors in the pseudogap regime 
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Experimental data obtained from tliermodynamic measurements in underdoped high temperature 
superconductors show unusual anomalies in the temperature dependence of the electronic specific 
heat both in the normal state and at the critical point associated to the superconducting phase 
transition. The observed deviations from the standard behavior are probably linked with the open- 
ing of a pseudogap in the energy spectrum of the single-particle excitations associated with the 
normal state. Based on a phenomenological description of the pseudogap phase we perform ana- 
lytical and numerical calculations for the temperature dependence of the specific heat for both the 
superconducting and normal state. The reduced specific heat jump at the transition point can be 
explained by a modified electronic single particle contribution to the specific heat in the presence of 
the normal state pseudogap. The hump observed in the normal state specific heat can be explained 
by the electronic pair contribution associated with strong fluctuations of the order parameter in the 
critical region. The obtained theoretical results are discussed in connection with experimental data 
for cuprates. 

PACS numbers: 74.20.Fg, 74.72.-h, 74.25. Bt 
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I. INTRODUCTION 



One of the most controversial properties of high tem- 
perature superconductor materials (HTSC) is the pres- 
ence of a gap in their normal state single particle excita- 
tion spectrum."'^ Usually addressed as the pseudogap, this 
relatively new feature just added a new controversy on 
the long list of unusual properties of the normal phase of 
HTSC. The pseudogap phase is seen in the underdoped 
region of the HTSC phase diagram for temperature val- 
ues above the superconducting critical temperature Tc 
and below a characteristic temperature T*. The pres- 
ence of the pseudogap phase was experimentally proved 
by direct measurements of the single particle excita- 
tion spectra in angle resolved photoemission spectroscopy 
(ARPES)^'- and tunnelling experiments'*, but also in nu- 
clear magnetic resonance (NMR), specific heat, resis- 
tivity, infrared conductivity, and Raman spectroscopy 
experiments^ However, in spite of a large amount of ex- 
perimental data, there is still no general consensus on the 
nature of the pseudogap phase, especially regarding the 
doping dependence of the onset temperature T* around 
the optimal doping point. As a function of doping, in 
one possible scenario, T* merges with Tc in the over- 
doped region, whereas in a different scenario, T* falls 
from large values in the underdoped region to a zero value 
at a critical point, universally identify for all HTSC at 
Per = 0.19.^ From the theoretical point of view, these two 
scenarios involve different approaches. In the first one, 
the superconducting and pseudogap phases arc strongly 
connected, the pseudogap being associated with the for- 
mation of precursor Cooper pairs (Tc < T < T*), pairs 
which become coherent and condense at Tc, leading to 



the superconducting phaseAl A different theoretical ap- 
proach leading to the same conclusions with respect to 
the onset temperature T* consider the role of the pair 
fluctuations above the critical temperature Tc^i^^ In 
the second scenario, the key role is played by the pres- 
ence of an antiferromagnetic region in the phase diagram 
at low doping values, the pseudogap being a consequence 
of the direct interaction between the electrons and fluctu- 
ations of the antiferromagnetic order parameter ^iiiS Un- 
fortunately, despite this large theoretical effort, an agreed 
description of the normal state in HTSC (including the 
pseudogap phase) is still lacking. 

In a recent work, Moca and Janko;^ performed a de- 
tailed theoretical analysis of the electronic specific heat 
in HTSC. Starting from a phenomenological description 
of the pseudogap, it is argued that a correct description 
of the specific heat behavior in the normal state can be 
obtained only with the inclusion of the electron pair con- 
tribution. Such a scenario is sustained by the presence 
of strongly enhanced fluctuations of the order parameter 
in the critical region above the transition temperature 
in quasi- two dimensional systems such as HTSC. In this 
way the observed maximumMii^ in the coefficient of the 
electronic heat capacity, 7(T) = C/T, can be fully ex- 
plained. However, the anomalies related to the specific 
heat behavior are observed also in the specific heat co- 
efficient jump, AjiTc), at the transition point:^ In the 
overdoped region, where the presence of a pseudogap 
is still questionable, the specific heat jump remains al- 
most constant. As the doping decreases, around the op- 
timal doping point, where the pseudogap is supposed to 
open, the specific heat jump starts to decrease. Once the 
doping value is in the underdoped region, AjiTc) falls 



2 



sharply, the larger the pseudogap is, the smaller the spe- 
cific heat jump becomes. The goal of this paper is to 
calculate the temperature dependence of the electronic 
specific heat, both for the superconducting and normal 
phase of HTSC, and to analyze its jump at the critical 
point based on a modified BCS theory which includes the 
pseudogap effects. A similar approach of the electronic 
specific heat behavior was considered by Loram et alm^ 
in order to explain the anomalous properties induced by 
the pseudogap in YBa2Cu306+a;- 

The paper is organized as follows: In Section II we 
discuss a phenomenological theoretical approach based 
on the validity of Gorkov's equations and we obtain the 
modified gap equation starting from a normal state char- 
acterized by the presence of a pseudogap, Eg. In Section 
III we discuss the specific heat behavior based on the pro- 
posed model. Analytical results are obtained for the spe- 
cific heat jump at the critical point using Pauli's theorem. 
Numerical results for the specific heat coefficient, 7(T), 
are presented based on an analytical expression of the 
free energy as function of temperature, below and above 
the critical temperature Tc. A comparison between an- 
alytical and numerical results for the specific heat jump 
at the transition point is also presented. Section IV gives 
our conclusions. 



II. THEORETICAL MODEL 

In the following we will consider a simple model based 
on the Gorkov's equations formalism in which the nor- 
mal state Green's functions will include the presence of 
the pseudogap. A similar analysis was used by different 
authorsiSiiLi^ in order to study the effect of the pseudo- 
gap phase on different properties of the superconducting 
state. Our theoretical approach is based on the assump- 
tion that in the pseudogap phase the self-energy correc- 
tions to the free electronic Green's function are given by: 



-£;2(k)Go(-k,-ic^„) 



(1) 



where Go(— k, — iw„) represents the free electron Green's 
function, £'g(k) the pseudogap, and a;„ = (2n -\- l)7rT 
is the usual fermionic Matsubara frequency. This phe- 
nomenological form of the self-energy was already used 
to explain the form of the spectral function, A{kp,u;), 
observed in ARPES experiments,'''' A similar behavior of 
the electronic self-energy was reported late in the seven- 
ties by Schmidt as a direct consequence of electron-pair 
fluctuation interactions in the critical region around the 
transition temperature. However, it is not our goal to un- 
derstand the origin of this phenomenological self-energy, 
but to use it in order to extract different properties of 
the superconducting state in cuprates. The normal state 
Green's function can be obtained with the aid of Dyson's 
equation as: 



G(k, iujj, 



(2) 



where El=^l + £;2(k), = (1 + ek/i?k)/2, and vl = 
(1 — Ck/£'k)/2 (Ck denotes the electron energy measured 
from the Fermi level). 

In terms of Green's function formalism the standard 
BCS theory is recovered by the use of the Gorkov equa- 
tions: 

A*{k)g(k,iLu^) - G-\-k,-iu;„)THk,iLu,,) = .(3) 

g{k, iujn) and J-{k, iujn) represent the normal and anoma- 
lous Green's functions in the superconducting state. The 
superconducting order parameter, A(k), is defined in the 
usual way in terms of the anomalous Green's function, 
T(k, icun), as: 

A(k) = -TY^ I J^Vik,p)Pip,^u;^) . (4) 

The interaction term, V(k, p), supposed to be attractive, 
is responsible for the formation of the Cooper pairs. The 
anomalous superconducting state Green's function can 
be easily obtained from Gorkov's equations and using 
Eq. Q the standard BCS gap equation is recovered. 

Our theoretical model assumes the validity of Gorkov's 
equations for the case of HTSC, where the free electron 
Green's function, Go(k, ia;„), is replaced by the more 
general Green's function given by Eq. In this way 

the effects of the pseudogap on the superconducting gap 
equation are considered. A simple calculation leads to 
the following general gap equation: 



A(k) = -r^ 



dp 

A(p) [0 



V{k,p) 



|A(p)P - Q] - _ ^2 _ ^2(p)] ' 



(5) 



The sum over the Matsubara frequencies can be per- 
formed analytically and the gap equation becomes: 



1 



dp 



A^{T) 



A4(r) + 4A2(r)i;2 



tanh 



B2(r) 



e^ + B^{TW{p) 



: tanh 



e2 + A2(r)^2(p) 
2r 



Q, + B^{TW{p) 

2r 



(6) 



where 



A^{T) = El - 
B^{T) = E^ 



A2(T) - 
A2(T) 



A4(r) + 4A2(r)£;2 



A4(r) -f4A2(r)£;2 



(7) 
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FIG. 1: (a) Numerical results for the temperature depen- 
dence of the superconducting gap function for different values 
of the ratio -E9/Ao(0) (Ao(0) represents the superconducting 
gap at r = in the absence of the pseudogap). (b) The zero 
temperature normalized superconducting gap, A(0)/Ao(0), 
as function of the ratio _Eg/Ao(0) obtained from both the 
numerical (full line) and analytical (dashed line) calculations. 



Unfortunately, the exact analytical solution for the gap 
equation ^ at any temperature in the superconduct- 
ing state (0 < T < Tc) is difficult to obtain. Instead 
we performed a numerical calculation for the gap tem- 
perature dependence considering different values of the 
normal state pseudogap. As a general observation we 
note that the superconducting state is suppressed by the 
presence of the normal state pseudogap. The value of 
the superconducting gap at T = K, A(0), transition 
temperature, Tj,, and superconducting gap, A(r), de- 
crease as the pseudogap increases. As an example of the 
numerical solution, in Fig. Q^a) we present the tempera- 
ture dependance of the superconducting order parameter, 
A(r), as function of temperature for different values of 
the pseudogap. Eg. 

However, some simplifications in Eq. © can be made 
in order to approximate some of the superconducting 
state properties, such as the value of the critical temper- 
ature, Tc, and the superconducting gap at T = 0, A(0). 
First of all the integration over the momenta is replaced 
by an energy integration using the corresponding density 
of states (DOS) in the normal state;^^i^^ 

where Nq denotes the DOS of a two dimensional (2D) 
free electron gas. This form of the density of states 
resembles the one seen experimentally in HTSC and is 
different from the constant density of states used in 
Ref. B We assume also the following energy scale. 
Eg < A(0) < W/2, with W being the bandwidth. This 
choice of the energy scale is valid in the doping region 
around the optimal doping point.- Based on this assump- 
tion one findsii^ 



ip{p) in the gap equation © is a factor associated with 
the general symmetry properties of the superconducting 
gap, pseudogap, and interaction potential. Experimen- 
tal data from ARPES and tunnelling experiments show 
that both the superconducting gap and pseudogap have 
the same symmetry, which in the case of HTSC is con- 
sidered to be of d-wave type^i^ with i/'Ip) — cos(20p) 
{9p = arctan(pj,/pa;)). Implicitly, the symmetry of the 
interaction term is assumed to be of the same type, i.e, 
y(p,k) — VdV'(k)V'(p)- The pure BCS case is simply 
recovered in the Eg^O limit. Note at this point the dif- 
ferences between our general gap equation © and the 
one used by Loram et al}^ in their approach. First of 
all, Eq. © includes contributions beyond the ones con- 
sidered in Ref. [3 which can be seen as a particular 
limit of Eq. © for small values of the pseudogap, EgM- 
Secondly, we considered a d-wave symmetry of the order 
parameter, which is known to be more appropriate for 
the case of HTSC.i-^ 



Tr_ = T. 



cO 



1 Eg 21C(3) / E, 



^ 

4rco 



A(0) = Ao(0) 



1 - 



E„ 



Ao(0) 
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3 
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a_ 



E„ 



Ao(0) 



(9) 



with Tco and Ao(0) being the critical temperature and 
zero temperature order parameter obtained in standard 
d-wave BCS calculations.-^ As we can see, the effect of 
the pseudogap is to decrease both the critical temper- 
ature and the superconducting gap at T = (See Ref. 
18), in agreement with the numerical results. Fig. ^b) 
compares the analytical result of Eq. © for the pseu- 
dogap dependence of the T = K order parameter with 
the exact numerical result obtained by solving the gap 
equation The slopes of the curves are slightly differ- 
ent, but such a behavior is expected as for the analytical 
results a number of simplifications were made. 
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III. SPECIFIC HEAT 

One of the most systematic experimental studies of the 
thermodynamic properties of the HTSC are the specific 
heat measurementsi^ii^. AnomaUes in the specific heat 
behavior are reported especially in the normal state, the 
coefficient of the electronic specific heat capacity, 7(T), in 
underdoped samples being no longer constant as function 
of temperature as we expect from the standard Fermi 
liquid theory. The broad maximum observed in 7(T) 
at a specific temperature Tm ^ T* is associated with 
the onset of the pseudogap. Moreover, the specific heat 
jump at the transition point is a function of doping, a 
lower value than the standard one predicted by the BCS 
theory being observed in the underdoped region. 

In this section we calculate the temperature depen- 
dence of the specific heat below and above the critical 
temperature, extract the specific heat jump at the critical 
point and compare our analytical and numerical results 
with the available experimental data. 



The electronic contribution to the specific heat can be 
calculated based on the thermodynamic potential as 
C = -Td'^Q/d'^T. One can see that the knowledge of 
the specific heat jump at the critical point requires the 
knowledge of the temperature dependence of the super- 
conducting gap in the critical region, below the phase 
transition point. A laborious, but straightforward calcu- 
lation based on Eq. gives 



A2(T) 



32^ 
21C(3) 



T 



1 - 3.14- 



Ao(0) 



7.6' 



Ao(0) 



(13) 



a result which inserted in Eq. H12|) gives for the specific 
heat jump at the transition point the following value 



A. Pauli's theorem 

According to the Pauli's theorem the difference be- 
tween the superconducting and normal state thermody- 
namic potential at the transition point can be calculated 
based on the following formula:^"' 



dA' (A') 



d{l/Vd) 
dA' 



(10) 



where 0s and 0„ are the thermodynamic potentials in 
the superconducting and normal state respectively, and 
V is the sample volume. In the critical region the di- 
rect dependence between the inverse of the interacting 
potential, 1/Vd, and the superconducting gap, A(T), can 
be extracted from Eq. Q considering that close to the 
transition point the value of the order parameter is small. 
Accordingly, one finds that 



1 



-T 



dp 



^p^p) (z0„)'-e 



-A^{T)Ty I 



(11) 



a result which together with Eq. H1U|) leads to 



&s - 0n 
V 

A4(r) 



dp 



(2^)2 



(12) 



1-0.53 



1.46 



( y 



21C(3) 



1 4- 2 62—^2 2 211 



(14) 



It is clear from Eq. (|14|) that the presence of the pseu- 
dogap in the single particle excitation spectrum in the 
normal state is responsible for the observed suppression 
of the specific heat jump at the critical point. Such an 
effect seems to be universal, as it is observed in the un- 
derdoped region of most of the HTSC materials phase 
diagrami^ii^'i^ 



B. Numerical results 

The electronic single particle contribution to the spe- 
cific heat has a simple linear T-dependence above Tc and 
an exponential temperature dependence below the crit- 
ical temperature. A detailed numerical analysis of the 
temperature dependence of the specific heat coefficient 
in underdoped HTSC was done in Refill. The main 
conclusion of this study was that the electronic single 
particle contribution to the specific heat is not enough to 
correctly describe the experimental data in the normal 
state, the inclusion of the electronic pair contribution 
being required in order to understand these data. The 
electronic pair contribution, associated to strong fluctu- 
ations of the order parameter in the critical region above 
the phase transition point, induces a hump in the normal 
state specific heat at a temperature of the order of T*, 
and has to be considered to understand the general be- 
havior of the specific heat in the normal state of HTSCii^ 
However, no numerical results were presented for the spe- 
cific heat jump at the transition point. 

In the following we will present numerical result for the 
temperature dependence of the specific heat coefficient, 
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7(T), in both the superconducting and normal state. Our 
analysis will include both electronic single particle and 
electronic pair contributions. From the numerical point 
of view, a more reliable way to extract the specific heat 
temperature dependence is to compute C as the first 
derivative of the electronic free energy, C — dE/dT. In 
general the electronic free energy, E, can be calculated 
based on the following relation:^^ 



E = T Y](ia;„ + ^k)5(k, ia;„) 



k,n 



Below the transition temperature the normal Green's 
function ^/(k, can be obtained from Gorkov equa- 
tions, while for temperatures above the critical tempera- 
ture the normal state Green's function, G(k, icjn), given 
by Eq ^ should be used instead. Two steps are per- 
formed in order to extract the temperature dependence 
of the specific heat coefficient. First, we will sum over the 
Matsubara frequencies and, secondly, we will numerically 
integrate over the momentum space, k. The general re- 
sult can be written as: 



E = Y,S^ 



(16) 



where 5k has different values in the superconducting and 
normal states. For the superconducting state, at tem- 
perature smaller the critical temperature, Tc, after the 
summation over the Matsubara frequencies is performed, 
can be calculated from the following relation 



5k = 



1 



i^2(k)./A(r)4 + 4A2(T)R2 



2^k^^(k) (A2(r) - El) VeITA(rmk) 



tanh 



A(r)2Vj2(k) 



2T 



2ek^'2(k) (52 (T) - El) .M^^mrWW) 



^(^)^^^«:tanhV^l 



B(r)2v,2(k) 



2T 



(17) 



whereas in for the normal state, at temperatures larger 
than the critical temperature, Tc, one has: 



5k 



1 



e2 + R2^2(k) 



2fk^/e^ + i?|^2(k) 



(2e^ + £;2V'2(k))tanh 



e2+i?2^2(k) 



2T 



,(18) 



The results obtained for the temperature dependence 
of the total energy, E, obtained by integrating Eqs. H17|) 
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Eg/Ao(0) 

FIG. 2: (a) Temperature dependence of the electronic sin- 
gle particle contribution to the specific heat coefficient, 7(r), 
for different values of the ratio iJg/Ao(0). (b) The nor- 
malized specific heat coefficient jump at the critical point, 
A7(rc)/A7o (Tc), (A70 (Tc) represents the value of the spe- 
cific heat coefiicient jump at the transition point in the ab- 
sence of the pseudogap) as function of the ratio i5g/Ao(0) 
obtained from both the numerical (full line) and analytical 
(dashed line) calculations. 



and (|18|l . match perfectly at the critical point but the 
slopes of the curves are different, indicating a disconti- 
nuity of the specific heat at the transition point. The 
specific heat coefficient can be calculated in a straight- 
forward manner by taking the derivative as a function 
of temperature. The results are presented in Fig. IHa). 
Two important features are observed: first, as the ana- 
lytical results predicted, the specific heat coefficient jump 
at the critical point decreases as the value of the pseudo- 
gap increases. Secondly, one can see that in the normal 
state, above the transition critical temperature, the pseu- 
dogap presence has no effect on the specific heat coeffi- 
cient, a constant value being obtained as the electronic 
single particle contribution is considered. Therefore, to 
explain the broad hump observed experimentally in the 
normal state specific heat, different mechanisms should 
be considered. Fig|2Ib) presents the normalized specific 
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- Electronic contribution 

- - Pair contribution 

- ■ Total contribution 



0.0 



0.5 



1.0 



1.5 



T/Tr 



2.0 2.5 



FIG. 3: Different contributions to the specific heat coefficient 
as function of temperature in the presence of the pseudogap 
(_Bg/Ao(0) — 0.4). The total electronic contribution (dotted 
line) to the specific heat coefficient is obtained as a sum of the 
electronic single particle contribution (full line) and electronic 
pair contribution (dashed line) . A simil ar qualitative behavior 
was reported for HTSC in Refs. llJl,! 



heat coefRcient jump at the transition point as function of 
the ratio £'g/Ao(0), as extracted from numerical and an- 
alytical calculations. As expected, there are slight differ- 
ences between the slopes of the numerical and analytical 
curves, a feature which is due to successive approxima- 
tions used in the analytical calculation. 

As we already mentioned different contributions to 
the specific heat coefficient, 7(T), have to be consid- 
ered for a better understanding of the HTSC specific 
heat data at high temperatures corresponding to the 
normal state. ^"^'^^ One possible mechanism was consid- 
ered by Moca and Jankoji^ and is related to the strong 
nature of the order parameter fluctuations in the criti- 
cal region, leading to an important electronic pair con- 
tribution to the specific heat coefficient in the normal 
state. This contribution has been evaluated previousl}*i^ 
and we quote here the result for completeness, C^""" ~ 
(r*/T)2exp(-2T*/T), where T* = ^nEj/g and g de- 
scribes the quasiparticle interaction. Fig ^ presents the 



temperature dependance of the specific heat coefficient 
in the presence of a normal state pseudogap {Eg/ A{0)). 
Electronic single particle and electronic pair contribu- 
tions are summed to obtained the total electronic specific 
heat coefficient. The electronic single particle contribu- 
tion is important for the specific heat coefficient jump at 
the transition point, whereas the electronic pair contri- 
bution can explain the broad maximum observed in the 
specific heat coefficient at high temperatures in the nor- 
mal state. Note that no specific heat coefficient jump at 
the transition point is associated to the electronic pair 
contribution. A similar qualitative behavior of the to- 
tal electronic specific heat coefficient was reported for 
HTSC^Jiii^ 



IV. CONCLUSIONS 

The nature of the pseudogap is still an open question 
and its explanation is beyond the scope of this paper. 
Using a simple phenomenological model we have inves- 
tigated the temperature dependence of the specific heat 
coefficient both below and above the critical point for 
HTSC characterized by the presence of a pseudogap in 
the excitation spectrum of the normal phase. The single 
particle contribution is the most important contribution 
to the electronic specific heat, but this is not sufficient to 
explain the hump that develops at high temperatures in 
the normal state. 

An important result that emerges from our calcula- 
tions is the behavior of the specific heat jump at the 
critical point. Experimentally, in the underdoped region, 
where the pseudogap energy is large, the jump is very 
small"" in contrast with the overdoped region where the 
pseudogap is practically absent and a BCS-like behavior 
emerges with a large jump at the critical point. We have 
calculated the electronic single particle contribution to 
the specific heat jump at the critical point in the pres- 
ence of the pseudogap and a similar dependence to the 
one experimentally observed was obtained. Increasing 
the pseudogap the specific heat jump at Tc starts to de- 
crease. This feature can be understand in terms involving 
the suppression of the superconducting state in the pres- 
ence of the pseudogap, as both the superconducting gap 
and critical temperature are smaller due to a loss of states 
at the Fermi level. The specific heat anomaly presented 
in HTSC at high temperatures (T > Tc) can be suc- 
cessfully understood in terms of two particle (electronic 
pair) contributions to the specific heat.«13, Note that the 
pair contribution to the specific heat does not affect the 
specific heat jump at the transition point. 

Our analysis is based both on analytical and numerical 
calculations. The restrictions imposed by the approxima- 
tions used in the analytical calculation lead to differences 
between the analytical and numerical results. Basically, 
the analytical results are reasonable around the optimal 
doping point and in the overdoped region, where the 
value of the observed pseudogap is small. With respect 
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to the phase diagram one should mention that the model 
we used is valid only in the second scenario discussed 
in the introduction, as the pseudogap presence leads to 
a reduction of the normal-superconducting phase tran- 
sition temperature, meaning that the pseudogap should 



completely disappear around the optimal doping point. 
However, it is generally accepted that in the overdoped 
region the standard Fermi liquid and BCS theories are 
still valid, any inclusions of a pseudogap in the descrip- 
tion of the system properties being inappropriate. 
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